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The avalanche region of an IMPATT (I M Pact ionization Avalanche 
Transit Time) diode under small signal conditions is characterized by the 
fraction of the total alternating current that is carried by holes and electrons 
in their respective drift spaces and by a residual impedance. The current 
fractions are roughly in phase with the total current below, and nearly 180° 
out of phase above a resonance frequency that is proportional to the square 
root of direct current density. 

This paper extends the calculations of Gilden and Hines for the current 
fraction to include phase shifts in the avalanche region so that extended 
avalanche regions can be considered. Realistic values (a ^ /3 in Si) for the 
ionization coefficients are used. Results of detailed numerical calculations 
for the current fractions as a function of frequency and direct current den- 
sity are presented. 

For typical frequencies and current densities, the residual impedance is 
negligible and hole and electron current fractions are equal. The avalanche 
region at a given frequency and current density is then characterized by 
one complex number and the admittance of a diode containing the avalanche 
region and adjacent drift regions is easily calculated. Plots showing the 
admittance as a function of frequency and current density for typical struc- 
tures are given. 

It is found that an optimal exponential growth rate of oscillations is ob- 
tained when the current density is such that the resonance frequency is about 
equal to one half the reciprocal transit time through the longest drift region. 
If the assumption is made that conditions giving the largest small-signal 
exponential growth rate give the best large-signal performance, an optimum 
Read-diode design is obtained for which the avalanche region width is a 
substantial fraction (tt $) of the drift region width. 

I. INTRODUCTION 

This paper considers the avalanche region of IMPATT* diodes, 1 - 2 - 3 
especially of Read'-type diodes in which the avalanche region is localized. 

* IMPact ionization Avalanche Transit Time. 

1797 
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Recent theoretical 4 and experimental 5 results for the impedance of Read 
diodes show considerable structure in the current and frequency depend- 
ence. This paper attempts to enhance the understanding of the small- 
signal negative resistance of IMPATT diodes by isolating the role of the 
avalanche region. Using the calculated parameters by which the ava- 
lanche region is characterized, the admittance for typical IMPATT 
diodes is calculated and exponential growth rate for oscillations is stud- 
ied. An application of the results to the design of Read diode oscillators 
is made. 

At a given angular frequency a and ac terminal current density 
t to t through the diode, the avalanche region can be characterized by 
three complex numbers: Z T , Fh , F e denned below. In Fig. 1, to the 
right of the avalanche region the small-signal particle current consists 
of a plane wave of holes. If we assume a constant drift velocity, then 
the magnitude is constant and the phase changes linearly with distance. 
Let Ah be the complex amplitude, extrapolated to a phase reference 
plane at x in the avalanche region. Then we define the hole current 
fraction Fh 



F h = Ah/itot ■ 



(1) 



Similarly, we define the amplitude of the electron current density, 
extrapolated to xo , to be A e and we define the electron current fraction 
F e 



F e = A e /i tot . 



(2) 



END 
REGION 
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Fig. 1 — Schematic of IMPATT diode. 
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In the limit as the frequency goes to zero, all the current through the 
device is particle current and F h and F e are of equal magnitude. The 
phase reference plane is chosen so that at low frequencies F h and F c 
are of equal phase also. The phase reference plane is approximately 
located at the point where the hole and electron currents are equal. 

Assume for a moment that the avalanche region were infinitely thin. 
From the fact that the total current equals the particle current plus 
displacement current we find that the ac field per total ac current, 
to the left of x is 



E 1 - i e /itot 1 - F e exp [iu(x - x )/v e ] 



(3) 



where i, is the electron current density and v e the electron velocity. 
The impedance, Z u between x h and x (see Fig. 1) is obtained by 
integration of (3). 



tO)C L L lUlT L J 

where t l is the transit time between .r and x L 



x Q — x L 

tl = 

V* 



(4) 



(5) 



and C L is the geometric capacitance per unit area, 

C L = _£— . (6) 

Xo — x L 

A corresponding expression obtains for Z 0R , the impedance between 

.To and x K . 

In deriving (4) it was assumed that the electron current is of the 

form 

i t = F e exp [ita(x — x )/v e ]- (7) 

This is true only outside of the avalanche region. The difference between 
the actual and asymptotic hole and electron current densities gives rise 
to an additional, residual impedance Z r . The total impedance, exclu- 
sive of end-region resistances, is then 

Z to t = Z r + Z L0 + Z 0l{ ■ (8) 

Typically, the residual impedance term is small in comparison to the 
other terms. Also, for narrow structures and at sufficiently low fre- 
quencies such that v/o) is large compared to the width of the avalanche 
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region, the hole and electron current fractions are equal in value and 
will be denoted by F without subscript. 

In Section II we consider, following Gilden and Hines, 6 the conse- 
quences of a simplified model of the avalanche zone, as embodied in 
Read's (13), (Ref. 1) and derive an expression for F from it 

_ 1 

F ~ 1 _(A 2 . (9) 



W 



Here the quantity f c = o> c /2ir, which we shall call the critical frequency, 
is real, independent of operating frequency a>/2ir and is proportional to 
the square root of the dc current density. 

Section III presents a treatment that is rigorous except for the 
idealizing assumption that hole and electron velocities are constant 
and equal in value. Numerical results for several structures are given 
in Section IV, and it will be shown that for frequencies / and direct 
current densities of interest F can be represented as 

F -(if-.A)(if-A)\=B)' (10) 

where J? is a zero on the real axis and is nearly independent of current 
density and A and A are complex conjugate poles that traverse a para- 
bolic path in the complex frequency plane as the current density is 
changed. 

In Section V a qualitative description of the change in negative 
resistance characteristics is given as a transition is made from p-i-n 
to p-n to "Read" diode. Admittance as function of frequency and cur- 
rent density is shown for two specific structures. Values for the expo- 
nential growth rate factor g (= — 1/2Q), maximized with respect to fre- 
quency, are shown as a function of current density. It is found that an 
optimal growth rate is obtained when the current density is such that 
the resonance frequency is about equal to one half the reciprocal transit 
time through the longest drift region. 

In Section VI we explore the consequences of the assumption that 
best large signal performance is obtained for the same conditions that 
yield largest small signal growth rate. This assumption in conjunction 
with a constraint expressing drift-region output limitations leads to 
Read-diode designs in which the width of the avalanche region is a 
substantial fraction (tt % ) of the drift region. 
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II. SIMPLIFIED MODEL 

The current fraction F can be deduced 6 from the simplified theory 
of Read. 1 His (13) for a thin avalanche region 

5* = /„(/<* -l) + 7. (ID 

states that the time derivative of the particle current density through 
the avalanche region, h , equals the product of I and a field dependent 
actor h(E) 

h(E) =Hjadx - l\ (12) 

if the saturation current I a can be neglected, as is usually the case. 
Under dc conditions h(E) = 0. If i p and e represent small signal ac 
components at frequency u of particle current density and electric 
field in the avalanche region, then 

iwip = h-—j-,e. (13) 

The total alternating current density i to t equals the sum of particle 
current density and displacement current density, 

itot = i P + *wee (14) 

or, with the definition 

F = i P /iiot , (15) 

.i- = -l (1 - F). (16) 

itot lue 

Division of (13) by i to t and substitution of e from (1(5) yields 

1 



F = 



1 -u " 
i\owit 



(17) 



with 



-u-atfit = y Jof|/e. (is) 

Equations (17) and (18) predict a pole in the impedance or a zero 
in the admittance at a critical frequency f CT \t that is proportional to the 
square root of the current density. Experimental measurements 7 and 
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numerical calculations 8 have shown that in p-n diodes and p-i-n diodes 
the admittance goes through a minimum, but not through zero, at 
frequencies approximately proportional to the square root of current 
density. This indicates that the simplified theory based on (11) is 
roughly valid but not accurate in detail. In the next Section a treatment 
is presented that does not make the assumptions inherent in (11) and 
that is, therefore, applicable to wide as well as narrow avalanche re- 
gions. 

III. IMPROVED MODEL 

In this section we give a more rigorous treatment of the avalanche 
zone. The present treatment is facilitated greatly by the idealizing 
assumption that hole and electron velocities are constant and equal in 
value. If this assumption is not made the problem is still tractable and 
a method of solution was obtained. 9 However, it is felt that a somewhat 
idealized treatment with the attendant reduction in complexity and, 
hopefully, gain in physical insight, is worthwhile. This is especially the 
case since experimental values for particle velocity vs electric field are 
available only for prebreakdown fields. Though electron velocities in 
silicon, the material of present greatest interest, can be extrapolated 
reasonably well into the breakdown region, considerable uncertainty 
prevails about hole velocities. However, since in silicon the electrons 
ionize much more strongly (roughly a factor of 10 more) than holes, it 
is felt that the results of calculations for the current fractions F are not 
affected significantly by the choice of hole velocity. Thus, the assump- 
tion that the hole velocity equals the electron velocity appears adequate. 
The value of 10 7 cm/sec is used for the velocity v. 

For the ionization coefficients a and (3 of holes and electrons as func- 
tion of electric field E the expressions 

a = 1.8 X 10 7 (cm -1 ) exp (-3.2 X 10 6 (V/cm)/E) 

P = 2.4 X 10 6 (cm -1 ) exp (-1.6 X 10 6 (V/cm)/E) 

are used. The numerical values refer to silicon at room temperature 
and are based on the work of Lee et al. 1 

Let //, and I e be the hole and electron current densities. Then the 
continuity equations state 

-^ 4 h = -h' + oJ h + pl e (19) 

v ot 

- d -L = U + ah + 01. (20) 

v dt 
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where primes denote spatial derivatives. Here we assume that holes 
move in the +x direction and electrons in the — x direction correspond- 
ing to a positive electric field as in an n-p structure. We introduce the 
total electric current associated with particle motion 

h = L + h (21) 

and the difference in currents 

h = I.- h. (22) 

Addition of (19) and (20) yields 

II Zj - I A ' + (0 + a)I z + (0 - «)Z A • (23) 

v dt 

From Poisson's equation 

h = -veE' + qvN D , (24) 

where E is the electric field and qN D the net impurity space charge. 
We shall now again consider small ac quantities at frequency w, de- 
noted by lower case symbols, superimposed on dc quantities and we 
denote by i tot the total ac current. From the continuity of total current 
we obtain 

*'tot = 4 + *'wee. (25) 

We define 

k = iu/v (26) 

and normalize the ac electric field to the total alternating current 



2 = — e. 

*tot 



(27) 



Then (25 ) can be written 



£ = 1 - kz. (28) 

^tot 

We introduce the derivative with respect to field of the last two terms 
of (23) 

1 



H = 

ve 



/.£(/^ «>+'*£<*-«>]. 



(29; 



Then the small signal ac version of (23) is 

-z" + (a + - k) (1 - kz) - - a)z + Hz - (30) 
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or 



[D- - Jc + (/3 + a)k + 03 - a)D - H)z = a + B - k, (31) 

where D is the spatial derivative operator, and where a, 0, and H are 
evaluated for the average, or dc, field as a function of distance. 

Equation (31) is a second-order differential equation for the complex 
quantity z. We postulate that to the left of the avalanche zone all 
particle current is electron current, i.e., 

is = plane wave moving to left (32) 

and that to the right of the avalanche zone all particle current is hole 
current 

ix = plane wave moving to right. (33 ) 

We do not specify that the magnitude of the ac electron current on the 
left equals that of the hole current on the right. Since we differ in this 
point from previous treatments, some discussion may be warranted. 

Holes and electrons are generated in pairs. One might, therefore, 
be led to conclude that the magnitudes of hole and electron particle 
currents must be equal. This is true for the dc or average currents but 
not generally so for the ac currents. Consider that a periodic generation 
rate g(x) is given. We allow g to be complex to represent the variation 
in phase with distance. The continuity equations (19) and (20) for 
the small signal case may be written 

(k + D)i h = g (34) 

(k - D)i. = g. (35) 

With g considered given, these are first-order differential equations for 
ih and i„ . With the boundary conditions (32 ) and (33 ) the solutions are 

ih{x) = I g(s) exp [— k(x — s)] ds (36) 

Jx L 

i„(x) = / g(s) exp [+k(x - s)] ds. (37) 

Since k is purely imaginary and g is complex, it is seen that ih(x R ) and 
i e {x L ) need not be the same. For narrow structures and at low fre- 
quencies, the difference in magnitude is negligible. This, however, is 
the result of the calculations and not imposed as a constraint. 

We specify, then, that outside the avalanche region the particle 
current consists of a plane wave moving to the left or right. 
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¥- - A R exp (-kx) right (38) 

*tot 

¥- = A L exp (+fcc). left (39) 

«tot 

The amplitudes A R and A L are as yet undetermined. By (28) the 
electric field outside the avalanche region is 

z _ \ [1 - An exp (-fee)] right (40) 

k 

z = j [1 - A L exp (+kx)] left (41) 

and we obtain the boundary conditions 

z' = 1 - kz right (42) 

z ' = _ (i _ fc 2 ) left (43) 

in accordance with (22) and (24). 

These boundary conditions specify completely the solution of (31). 
If the solution for z is carried over the entire depletion region, then the 
integral over z gives the total impedance, exclusive of resistive losses 
in the end regions. However, in order to isolate the role of the ava- 
lanche region and to bring the impedance into the form (8) we proceed 
as follows: If a solution for z has been obtained at some sufficiently 
low frequency, then the magnitudes of A R and A L in (38) and (39) 
are equal. A value x exists at which the right hand side of (38) equals 
the right hand side of (39); i.e., at which the asymptotic particle cur- 
rents, extrapolated into the avalanche zone, are equal in phase as well 
as magnitude. We call x the phase reference point and define the 
extrapolated particle currents at this point as the current fractions 
F e and F h ; 

F h = A H exp (-fca>) (44) 

F e = A L exp (+fcr ). ( 4 5) 

We define the asymptotic particle current i a , 

(F h exp [— k(x — .To)] x ^ x 
iji^ = (46) 

[F e exp [+k(x - .To)] x < .To 

and the residual particle current i T as the difference between the actual 
particle current ?" to t(l - kz) and the asymptotic particle current 
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^- = 1 - kz - ^- . 



«tot 



*tot 



(47) 



Note that by construction i r vanishes outside the avalanche zone and 
is thus independent of the dimensions of the drift zones, as are F e and 
Fh . We can now solve (47 ) for z, 



1 

Z = k 



*tot 



(X exp [— k(x — x )] 
[F 9 exp [+k(x - x )]J 



X ^ .To 



X < Xq. 



(48) 



The integral of z, multiplied by l/ve = k/ (uoe) is the total impedance 
Z = Z r + X ^^ \\ - F e ~ exp (iu>TL) ] 



litiE 



= Z r -\- Z L o + z 0R 
with 



+ 



Xr — Xq 



[i-fti 



— exp (itoTjj) 



>„=-± (±dx. 

l(i)£ J lint 



_ (49) 



(50) 



(51) 



We shall use the symbol F without subscript to refer to either F e or 
Fh if their difference is negligible or the distinction unimportant. 



IV. NUMERICAL RESULTS 

This section contains numerical results for the following structures in 
silicon : 

D3 error function complement p-diffusion from a surface con- 
centration of 10 20 cm -3 into a 3 X 10 16 cm -3 n-type substrate 
with a junction depth of 3 microns. 

D12 similar to D3, but with a 12-micron junction depth. 

L22 a linearly graded junction with a concentration gradient of 
10 22 impurities/cm . 

C0.1 Constant field avalanche zone (field = 5.71 X 10 5 volts/cm) 
of 0.1-micron width, surrounded by regions of sufficiently low 
field that negligible avalanching takes place there. 

CI Similar to C0.1, but with 1-micron width; field = 3.57 X 10 B 
volts/ cm. 

C5 Similar to C0.1, but with 5-micron width; field = 2.79 X 10 6 
volts/cm. 
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CIO Similar to C0.1, but with 10-micron width; field = 2.52 X 10 6 
volts/cm. 

For each desired direct current density h , self-consistent values of 
dc electric field E, difference current density h , ionization coefficients 
a and /3, and derivative quantity H, (29), were computed as described 
in Appendix A. Using these quantities, the differential equation (31), 
subject to boundary conditions (42) and (43), was solved numerically 
for ^-values corresponding to frequencies of interest. 

Figs. 2 and 3 show the results for the quantity 1 — l/F = ( — dis- 
placement current/particle current) for direct current densities of 100 
and 1000 amps/cm 2 . According to the simplified model of Section II, 
1 - l/F should have a real part varying as the square of frequency 
and a vanishing imaginary part. As Figs. 2 and 3 show, the square 
law is obeyed quite well at low frequencies and current densities. Only 
for the wide structures and/or at high frequencies and current densities 
does F deviate from the square law and do F c and F h deviate from each 
other appreciably. The frequency / crit at which (1 — \/F) reB .\ is unity 



10 

8 

b 















'// / 




















i 




















'/&/ 
















A -V 


y/? 










iooA/cm 2 






Vl. 


/ 














/ i 


' / 


// 
















/ / 


// 


// 
















/ / 


// 


'/ 


















/A 


1 














// 


/ If 


//_ 
















/ / 
















/ 


y 


i / // 
















/ 




/ // 
















// 


/ 


Jr_ 














/ 


/, 


z /> 


7 















I 10 

FREQUENCY IN GHZ 



Fig. 2 — Real part of 1 - l/F as function of frequency for various structures. 
Current density 100 amps/cm 2 . 
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Fig. 3 — Real part of 1 — 1/F as function of frequency for various structures. 
Current density 1000 amps/cm 2 . 

is in good agreement with values obtained from the simple theory of 
Section II. 

The absolute value of the imaginary part of 1 — 1/F at a current 
density of 100 amp/cm 2 is plotted in Fig. 4 as a function of frequency / 
for the various structures. The current dependence of the diffused diode 
D3 is shown in Fig. 5; that for the other structures is similar. The 
vertical line on each curve to the right of the maximum indicates the 
frequency at which the real part goes to zero. At frequencies 15 to 25 
percent higher, (1 — l/F) imaK goes through zero, being positive for 
lower and negative for higher frequencies. Asymptotically, at low 
frequencies (1 — 1/F)in iag is proportional to +/ and at higher frequen- 
cies to — f 3 . The similarity of the curves is striking and suggests that a 
fit with a few parameters ought to be possible. The most meaningful 
expansion, and one working well here, is to consider F as an analytical 
function of a complex frequency variable and to expand it in terms of 
the poles and zeroes nearest the frequency domain of interest. At cur- 
rent densities of interest and for frequencies / from zero to a few times 
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/ ( , rit , the numerical results can be represented well by two poles and a 
zero; 

(if - B) AA 



F = 



(52) 



(if- A)(if - A) -B' 

where ~ denotes the complex conjugate. A is the complex resonance 
frequency 

A=f + if r . (53) 

We denote 

B= -f c . 

We utilize the following inequalities which hold for current densities of 
interest: 



/.«/r 
ft « fc 



(54) 
(55) 
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Fig. 4 — Imaginary part of 1 — 1/F as function of frequency for various struc- 
fcures. Current density 100 amps/cm*. The vertical bars indicate /crit. 
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Fig. 5 — Imaginary part of 1 — 1/F as a function of current density for diode 
D3 at 100 and 1000 amps/cm 3 . 



and we consider frequencies small compared to f c . Then the reciprocal 
of (52) can be written 



1 



= 1 



J-*C. + $-j&)- (56) 



7. 

Equation (56) has the frequency dependence embodied in Figs. 2 
through 5, i.e., the real part of 1 — 1/F changes as frequency squared 
and the imaginary part is the difference of a linear and cubic term in 
frequency. The low-frequency asymptote of the imaginary part is 



(1 - 1/F) ,„„=///., 



(57) 



with 



/. = 



l + 2fjc/fr 



(58) 



The frequency, / cr it , at which the real part of 1/F vanishes equals 



fr- 
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font = /, (59) 

and the frequency/, at which the imaginary part of \/F vanishes is 

/, = Vfr* + 2fafc ■ (60) 

Now let us consider the dependence of F on the direct current density 
/. Three relations hold approximately for the curves of Figs. 2 through 

5. 

The critical frequency / cr it is proportional 

to the square root of direct current density. (61) 
The low-frequency, linear asymptote of (1 — l/F) imw 

is independent of current density. (02) 
The ratio of the frequencies / crit and /,- at which 

the real and imaginary parts of \/F cross zero 

is independent of current density. (63) 
We show now that these relations lead to the following current 
density dependences : 

u - VI < 65 ) 

f c independent of I. (66) 

Equation (65) follows directly from (59) and (61). Taking the ratio of 
(60) and (59) and using (63) we find that 

fafc/fr 2 is independent of /. (67 ) 

Use of (62), (57), (58), and (67) yields (66). Finally, from (65), 
(66), and (67) we obtain (64). Thus, as long as the inequalities (54) 
and (55) and relations (61) through (63) hold, the current density 
dependence of F can be expressed as follows (see Fig. 6): With increas- 
ing current density the poles A and A of F move through parabolas, starting 
at the origin, while the zero at -f c is independent of current density. Where 
F e and F h deviate from each other, the representation (52) is still usable. 
The poles A and A are the same, but F e and F h have separate zeros at 

— f et and -fch ■ 

The above discussion shows how the poles and zero can be obtained 
from calculated results of F as a function of frequency and current den- 
sity. In Appendix B a more direct way of obtaining the poles as solutions 
of an eigen -value problem is presented. The results are shown in Fig. 6, 
where the poles are plotted in the complex frequency plane with cur- 
rent density varying along the curves. The zeroes are listed in the insert. 
Only positive imaginary values are shown. Symmetry exists about the 
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FREQUENCY, f a IN GHZ 



Fig. 6 — Locus of poles of F as function of current density in the complex 
frequency plane. Note unequal scales of real and imaginary axes. 

real axis. Note that the scales of real and imaginary axes differ by a 
factor of 10. 



V. NEGATIVE RESISTANCE AND MAXIMUM GROWTH FACTOR 

In this section, we consider the interaction of the avalanche and drift 
regions in causing negative resistance. The most important aspect of 
the current fraction F is the rapid change of phase, with respect to total 
current, from near zero below, to more than 180° above, the resonance 
frequency. This behavior follows from (52) and is illustrated for struc- 
ture D3 in Fig. 7 where the phase p P of F is shown for current densities 
of 100 and 1000 A /cm 2 . At the higher current density the transition is 
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more gradual. Likewise, for wider structures the transition is more 
gradual than for narrower structures. 

Let us consider a diode in which the hole drift region is negligible 
and for which the residual impedance Z T is negligible so that the dom- 
inant term in the impedance is 

— exp (— iWl) 






IdiTL 



']■ 



(68) 



We want to explore under what conditions the small signal impedance 
has a negative real part. This will occur when the quantity 



„ 1 — eXp (-tMTl) m pg 



(69) 



has a phase angle between and 180°. <p a , the negative of the phase of 
G, is plotted in Fig. 8 as a function of frequency. We define as drift 
frequency, /drift , the frequency at which the phase (p a is 90°, 



/drift = 



2(z — x L ) 



2r L " 



(70) 



The singularity at twice the drift frequency occurs where G goes to 
zero and the phase is indeterminate. If the operating frequency is below 
1/tl we refer to the diode as operating in the ir mode. The mode- 
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Fig. 7 — Phase of F for diode D3 at 100 and 1000 amps/cm 2 . 
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Fig. 8 — Negative of drift delay phase. 



designation for higher frequencies is given in Fig. 8. The condition on 
the phases for negative resistance can now be stated 



< <p, - <p a < 180 c 



(71) 



and can easily be visualized if Figs. 7 and 8 are combined. This has been 
done in Fig. 9 where a drift frequency of 60 GHz appropriate for the 
depletion layer of diode D3 has been used. Starting with the 100 A/cm 2 
curve for <p F , we see that negative resistance prevails over a wide 
frequency range, as indicated by the shaded region. Near/ = 2f dri{t 
the magnitude of G is small, hence the negative resistance contribution 
of FG is small and is likely to be outweighed by parasitic end resistance. 
At frequencies beyond 2/d r if t negative resistance occurs in higher-order 
modes. At frequencies much higher than / crit the phase tp P depends on 
frequency in a complicated way and the present calculations lose ac- 
curacy. As the current density is increased, the <p P curve moves to the 
right with some softening of the transition. The lowest frequency at 
which negative resistance sets in increases also, but a wide range of 
negative resistance continues to exist. The behavior thus far described 
is typical of a single-diffused p-n junction diode. 

Next, consider a typical "Read" structure, i.e., one in which the 
drift space is much wider than the avalanche region. We consider the 
same avalanche zone as before but let the drift space be, say 8 microns 
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with a drift frequency of 6.25 GHz, as shown in Fig. 10. We now see 
that the frequency range for negative resistance in the it mode is much 
narrower; at 1000 A/cm 2 there is no r-mode negative resistance for this 

structure. 

Now let us consider structures with narrow drift regions. The extreme 
case for which the drift regions coincide with the avalanche region 
is represented by p-i-n diodes. 8 All three terms in (8) contribute, but 
the conclusions reached from a study of Z L0 alone are qualitatively 
correct. Negative resistance is obtained for a wide range of frequencies. 
For a consideration of the low-frequency behavior we expand F to 
first order in / 

(72) 



73) 







F w 1 + i f 7 




(cf. (57) and (58). 


With 








G& 


i _ *p = i _ 


2/,irirt 


(68) yields 










Rlo = 


• ' d- 


2/<l r ift\ 




4/ <lt „, C L V 


Thus, if 




/drift < Jj/» J 





(74) 



(75) 

negative resistance prevails from it mode frequencies down to dc. 
The case of the uniform-field avalanching plasma has been discussed 
by Misawa. 

For an interpretation of the frequency /« , consider the change in 
electric field 8E, per change in total current in the limit as the frequency 
goes to zero, as shown in Fig. 11, for diode D3 at 100 A/cm 2 . The space 
charge associated with the additional current causes 8E to be roughly 
hyperbolically shaped with the apex at the point where hole and electron 
currents are equal. The placement of the 8E curve with respect to zero 
is such that reduction of ionization in the center where 8E is negative is 
compensated by an increase in the adjacent regions. Outside the ava- 
lanche zone 8E changes linearly. At points A and B, separated by dis- 
tance d AB the asymptotes to 8E outside the avalanche zone cross zero. 
From a consideration of (28) through (32) it can be seen that l// f is 
half the transit time between points A and B, 
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Fig. 9 — Superposition of Figs. 8 and 9 appropriate for p-n junction case. 
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Fig. 10 — Superposition of Figs. 8 and 9 appropriate for "Read" diode case. 
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P-SIDE 




Fig. 11 — Incremental field at low frequencies for p-n diode. 



, _ cIab 



(76) 



If the avalanche zone is narrow compared to the drift region the integral 
over 8E is positive and so is the dc incremental resistance. 

For a p-i-n diode (CI at 100 A/cm 2 ) the incremental field 8E is as 
shown in Fig. 12. The asymmetry results from the inequality of hole 
and electron ionization coefficients. It is seen by inspection that the 
integral over 8K is negative and also that the transit time between A 
and B is less than tt/2 times the transit time through the depletion 
layer (this is the condition corresponding to (75) when both left and 
right drift regions are considered), hence the resistance is negative. 
It is to be noted, however, that if equal hole and ionization rates had 
been used in this calculation, zero incremental dc resistance would 
have been obtained. 

Up to this point we have considered the phase conditions for the oc- 
currence of negative resistance. Now let us consider the actual resistance 
or, more conveniently, the admittance Y. Fig. 13 shows plots of the 
imaginary vs real part of Y, with frequency varying along each curve, 
for diode D3, area = 10 -4 cm 2 , at various current densities. Fig. 14 
shows a corresponding plot for a Read structure Rl consisting of ava- 
lanche region CI (1 micron wide, constant field) and an adjacent 
electron drift space 9 microns wide. The diode area is also taken as 
10 -4 cm 2 . The admittances in Figs. 13 and 14 are the reciprocal total 
impedance, but differ insignificantly from what would have been ob- 
tained with neglect of the residual impedance. 

Now let us consider the quality factor Q or the related growth factor g, 
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Fig. 12 — Incremental field at low frequencies for p-i-n diode. 
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' - - k (77) 

which has the following meaning: If a diode having negative resistance 
is imbedded in a lossless circuit in which it oscillates at frequency / = 
co/2t, then the absolute values of voltage and current amplitudes vary 
with time as exp (gut). It is convenient to introduce the complex fre- 
quency variable s 

s = (g + i)f. (78) 

From (4), (52), and (53) it can be inferred that under open-circuit 
conditions the diode will oscillate at a value 

8 - A = /. + ifr 

and, therefore, the open-circuit growth factor is 

Qopen = Ja/Jr • 

According to the current-dependence of Fig. 6, g ope n is proportional to 
the square root of the current density. Note that under open-circuit 
conditions g depends on the avalanche structure only and is independent 
of the drift regions. This is not the case under more general circuit 
conditions. 

If the admittance Y is an analytic function of frequency, as is physi- 
cally reasonable and as we shall assume to be the case here, then the 
admittance is denned also for complex frequencies. The largest growth 
rate g (/) for reactive circuit imbedding at a given frequency / is that g 
for which 

m + i)f] + arte + m - ° (79) 

if Y is capacitive; or 

Y[(g + i)f] + 2ir(g + i)fC - (80) 

if Y is inductive. Here L or C are chosen so as to resonate the imaginary 
part of Y. The quantity g thus obtained is related to the quality factor 
Q, as conventionally denned, by (77). So far, g Ls defined at a given 
frequency / and current density. We define as g max the maximum of g 
(or maxima, where relative maxima exist) with respect to frequency 
at a given current density. Finally, we define as g opt the maximum of 
fifmax with respect to current density. We define by / max the frequency (ies) 
at which max occurs. Figs. 15 and 16 show g max and f max as a function 
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Fig. 15 — Maximum growth rate and open-circuit growth rate of diode D3. 

of current density for diode D3 and Figs. 17 and 18 for diode Rl. Also 
shown for comparison are f r and ^ opcn . Included in Figs. 16 and 18 are 
horizontal lines marked / dr ift representing the drift frequency of the 
longest (= electron) drift region. An important result contained in 
Figs. 15 to 18 is that g oP t is obtained for current densities at which 
the resonance frequency f T is about equal to the drift frequency and that 
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Fig. 16 — Frequency at which gr mo x occurs and j r for diode D3. 
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Fig. 17 — Maximum growth rate and open circuit growth rate for Read diode 



/max at this current density is some 20 to 30 percent above f T . For 
the diffused diode D3 the drift frequency is high (tt 60 GHz) and there- 
fore, g opt is obtained at a high current density, 3 X 10 4 amps/cm 2 . 
This current density is so high that a substantial widening of the ava- 
lanche region, with respect to its low-current-density configuration, 
has taken place. As a consequence, f c now is no longer independent of 
current density, but rather decreases with current density. Deviations 
from the low-current-density relations are also seen in the break- 
away from the line of slope \ for ^ oP en in Fig. 15. 

For the Read diode Rl the largest g op t , occurring in the tt mode, is 
obtained at 200 amps/cm 2 . The values of g max in the nr mode fall off 
rapidly as the frequency approaches (n + 1 )/drift • For a Read-type 
structure, i.e., one having a substantial drift space in which negligible 
avalanche multiplication takes place, g op t is near 0.13 and is only 
weakly dependent of the detailed structure of the avalanche region. 
On the other hand, for diodes in which the avalanche region occupies 
all of, or a substantial fraction of, the depletion region g opt is closely 
related to the value of ^ op en and is larger than 0.13, as is the case for 
D3. 



VI. LARGE SIGNAL DESIGN CONSIDERATIONS FOR READ DIODES 

In Section V we obtained the result that the optimum growth rate 
under small-signal conditions occurs when /drift ~ f r • This condition 
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provides a relation between current density, drift region width and 
avalanche region width (via/ r ). It is plausible to assume that the same 
relation should be approximately satisfied for best large-signal perform- 
ance. Preliminary large-signal results* have shown that a direct corre- 
spondence between the large-signal and small-signal properties does not 
exist and that, for example, large-signal self-sustained oscillation can be 
obtained for frequency-current-density combinations for which the 
small-signal resistance is positive. Nevertheless, as a point of departure 
for a large-signal design we choose a structure in which avalanche region, 
drift region, and current density are so related that g opt is obtained at 




6 8 10 4 
CURRENT DENSITY IN A/cm 2 

Fig. 18 — Frequency at which g max occurs and j r for diode Rl. 

the desired frequency under small signal conditions. This provides one 
constraint. 

Resonance frequency/, as a function of current density for the various 
structures is shown in Fig. 19. Suppose a 5-GHz operating frequency is 
wanted. We would need a resonance frequency about 25 percent lower 
or 4 GHz. Thus, structure CI at 120 amps/cm 2 , or C5 at 400 amps/cm 2 , 
etc. could be used. The drift region would be 12.5 microns wide. 

Another constraint between average current density and frequency is 
imposed by the output limitations of the drift space. The following 
discussion applies to Read diodes having a constant-field drift region 
only. As shown by Read, 1 under large-signal conditions the carrier cur- 
rent through the drift region is carried in the form of charge pulses. If 

* D. L. Scharfcttcr and H. K. Gummel, work in progress. 
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/ is the frequency of operation and I DC the average current density, 
then the charge per pulse is 



This pulse causes a change in electric field of 



AE = 



Idc 



(81) 



(82) 



Let Adrift be the field in the drift region at the onset of breakdown; 
we shall use a value of 1.5 X 10 V/cm which is reasonable for silicon. 
AjE7 must not be larger than approximately •£ £"drift • The factor of ^ 
takes into account that in the desired mode of operation the terminal 
voltage reaches its minimum while the charge moves through the drift 
region and that therefore an additional lowering of the drift field below 
Adrift occurs. Thus, we have the condition 



I DC < \£E dim f- 



(83) 



This constraint is shown in Fig. 19 by the line of slope 1. Only the region 
to the upper left is allowed. For a diode to have a large power capability 
the highest feasible current density should be used. Thus, a design cor- 
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Fig. 19 — Illustration of Read diode design considerations. 
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responding to the circle appears most promising. The diode would operate 
at a current density of 300 amps/cm 2 and have an avalanche region 4 
microns wide. Note that the avalanche region width would be about 
| of the drift region width. This is in contrast to Read's proposal that 
the avalanche region should be very narrow in comparison to the drift 
region. However, the design suggested here is based on a combination 
of small-signal dynamics and large-signal output limitations and may 
have to be modified when large-signal results become available. 

appendix A 

This appendix describes the steps by which self-consistent dc solu- 
tions for electric field E and difference current density h are obtained 
for a given terminal current density is . The electric field is the sum of 
a component E D (x), related to the space charge of the impurities and a 
component 

E c (x) = (&/h)z(x), 

where z' = h/h , related to the space charge of the mobile carriers. 

(t) Select a point xi in the avalanche region and assign a trial value 
#i for the electric field at this point. Conveniently, Xi may be the loca- 
tion of the metallurgical junction, but the choice is not critical. 

(it) Choose a trial function for the difference current h . A conveni- 
ent, though crude, choice is 

h = h , x £ X! (84) 

7a = -h , x > x x . (85) 

(Hi) Compute the electric field E 

E(x) = Ex + I jf [^ + qN D (a)~\ ds. (86) 

Extend the integral up to boundary points x h and x R to the left and 
right of .Tx , i.e., to points where the field reaches the value required to 
carry the current through the not swept-out semiconductor region : 

E(x L ) = Iz/W I N D (x L ) | ] (87) 

EM = /z/U, I NdM I ]. (88) 

Define 

E D (x) = Ei + q -[ N D (s)ds (89) 
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Use 

z (x) = f I A (s)/h ds (90) 

J xi 

as trial function for z. The beginning of the next step, (iv), is redundant, 
but is required for the iteration loop. 

(iv) Set E(x) = E D (x) + (ev/h)z (x), and /a = hz '. Using these 
values, evaluate a, fi, and H, (29). 

(v) Solve 

[D 2 + (/3 - a) D - H]z = a + 13 - Hz (91) 

for z with boundary conditions z' = — 1 at left boundary and z = + 1 
at right boundary. 

(vi) Make z the new trial solution z . 

(vii) Repeat steps (iv) through (vi) until z and z differ negligibly. 
If the final field at x\ differs appreciably from the trial field E\ the bound- 
aries are incorrectly defined, and steps (i) and following should be re- 
peated with the current value E(xi) as trial field Ei . This time use /a = 
z'Ts instead ol (84) and (85) as trial function for / in step (ii) . 

(viii) When the conditions of step (vii) are satisfied, i.e., when z tt z 
and E(xi) tt Ei , then E, I& = z'l? , a, /3, and H are the desired self- 
consistent dc quantities. 

APPENDIX B 

In this appendix we show how the complex resonance frequency 
A = f a + ifr at which the impedance goes to infinity can be obtained 
directly from (31), rather than by curve fitting of results obtained at a 
set of imaginary s-values s = + if, or real frequencies. At resonance 
the normalized field z goes to infinity and therefore the right-hand side 
of (31) is negligible compared to the left-hand side, i.e., we have the 
homogeneous equation 

[D 2 - k 2 + (a + 0)k + (a - /3)D - H]z = 0. (92) 

Solutions of (92) with boundary conditions (42) and (43) exist only for 
special values of k. The quantity k(k — a — 0) may be considered the 
eigenvalue of (92). Of the possible k values we select the complex 
conjugate pair of smallest absolute value. Let k be the k with positive 
imaginary part. Then 

vh = 2ir(f a + ifr) (93) 

is the desired resonance frequency. 
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Numerical values of k may be obtained as follows. The differential 
equation (92) is replaced by a set of appropriate difference equations 
for the values z n of z at a set of meshpoints. This set of equations may 
be written in matrix form 

M • Z = 0, (94) 

where Z is a vector having the z„ for its elements. The matrix M, being 
derived from a second-order differential equation, is a "tridiagonal 
matrix", i.e., it has nonvanishing elements only on the main diagonal 
and the two adjacent diagonals. Now k must be chosen so that the 
determinant of M vanishes. If fci and fc 2 are trial values and the cor- 
responding values of the determinant are Di and D 2 , then one may 
choose as next trial value 

S , _ _ Wfe-W (95) 

and thus iteratively approach the desired value of k • 

The frequency f c of the zero is chosen so that F has the proper phase 
in the neighborhood of the resonance frequency. Let F be the current 
fraction at the complex frequency Si = /« + ifi , where /] is near /, , 
then 

/c=/i^-/ Q . (96) 

' imag 
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